1 Method for Correcting the Effects of Interdetector Band Broadening. 

2 

3 Related patents 

4 The methods disclosed in this specification will have particular relevance to a variety of patents 

5 and co-pending applications by the assignee and inventor. They include the following: 
6 

7 Steven P. Trainoff and Philip J. Wyatt, Application 10/205,637 filed 24 July 2002, "Method for 

8 determining average solution properties of macromolecules by the injection method." 
9 

1 0 Philip J. Wyatt, U. S. Patent Number 6,41 1 ,383 25, issued 25 June 2002, "Method for measuring 

1 1 the 2 nd virial coefficient." 

12 Steven P. Trainoff, U. S. Patent Number 6,180,906, issued 30 January 2001, "Electrode design 

13 for electrical field flow fractionation." 
14 

15 Steven P. Trainoff, et al., U. S. Patent Number 6,128,080, issued 3 October 2000, "Extended 

16 range interferometric refractometer." 
17 

18 Gary R Janik, et al, U. S. Patent Number 5,900,152, issued 4 May 1999, "Apparatus to reduce 

1 9 inhomogeneities in optical flow cells." 
20 

2 1 Gary R Janik, et al, U. S. Patent Number 5.676,830, issued 1 4 October 1 997, "Method and 

22 apparatus for reducing band broadening in chromatographic detectors." 
23 
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1 Philip J. Wyatt, et al, U. S. Patent Number 5,530,540, issued 25 June 1 996, "Light scattering 

2 measurement cell for very small volumes." 
3 

4 David W. Shortt, U. S. Patent Number 5,528,366, issued 1 8 June 1 996, "Precision determination 

5 for molecular weights." 
6 

7 Philip J. Wyatt, U. S. Patent Number 5,305,07 1 , issued 1 9 April 1 994, "Differential 

8 refractometer." 
9 

10 Background 

1 1 The analysis of macromolecular species in solution by liquid chromatographic techniques is 

1 2 usually achieved by preparing a sample in an appropriate solvent and then inj ecting an aliquot 

1 3 thereof into a chromatograph. The chromatograph includes fractionation devices of various types 

14 that separate the sample as it passes through them. Once separated by such means, generally 

1 5 based on size, mass, or column affinity, the samples are subjected to analysis by means of light 

16 scattering, refractive index, UV, viscometric response, etc. For example, in order to determine 

1 7 the mass and size distributions of a particular sample whose separation is performed by size 

1 8 exclusion columns, the chromatograph would measure the sample sequentially by multiangle 

1 9 light scattering, MALS followed by a differential refractometer, dRI. The dRI determines the 

20 concentration while the MALS unit measures the excess Rayleigh ratio as a function of angle of 

2 1 each eluting fraction. For molecules very much smaller than the wavelength of the incident light, 

22 light scattering measurements at a single angle may often be sufficient. 
23 
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As each fraction passes through a detection device, it produces a signal often referred to as a 
"band" or "peak." Because of dispersion and mixing effects, these bands are broadened 
somewhat each time the sample passes through a different device. Consider a sample comprised 
of a low concentration aliquot of a monodisperse protein. In this event, the Rayleigh excess ratio 
is directly proportional to the molar mass and the concentration. The light scattering signal and 
the concentration signal should be of identical shape and would overlay perfectly were the two 
responses normalized to have the same areas. However, as the sample passes from the MALS 
detector and enters the dRI detector, it passes through intermediate regions and connections that 
contribute to the dispersion and mixing of the sample. Many sources of band broadening and 
some of their remedies are discussed in detail by Yau, et al. in their book "Modern size exclusion 
liquid chromatography" published by John Wiley & Sons in 1979. In the above example, the dRI 
signal will always appear somewhat broadened with respect to the MALS signals. 

The effect of broadening is shown clearly in Figure 1 which shows a chromatogram of a sample 
of bovine serum albumen, BSA, in an aqueous buffer along with the computed molar mass, 
uncorrected for the effects of broadening. The uncorrected molar mass is shown in trace 1. It is 
presented divided by the molar mass of the monomer, M w0 , so that the vertical axis for the 
monomer reads 1, the dimer reads 2, etc. Trace 2 corresponds to the 90° light scattering signal as 
a function of time. Trace 3 corresponds to the dRI signal. The aggregation states have been 
fractionated so that the rightmost peak region, from 16 minutes to 18 minutes, is data from the 
pure monomer, from 14 minutes to 16 minutes is from the dimer, etc. In each peak, broadening 
causes the concentration at the center of the peak to be suppressed and that in the "wings" to be 
enhanced. Combining this broadened dRI signal with the MALS signals to determine molar mass 
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1 at each elution time causes the molar mass determined near the peak's center to be systematically 

2 overestimated, while that in the wings to be underestimated. The situation reverses when the 

3 concentration detector is upstream of the MALS detector. This gives rise to the non-constant 

4 molar mass data within each peak region. If there were no band broadening, the data would 

5 consist of a series of plateaus at 1 for the monomer, 2 for the dimer, etc. As the fractionated 

6 sample becomes more polydisperse, the effects of band broadening become less apparent (but are 

7 still present). For a well fractionated sample comprised of discrete aggregate states, the problem 

8 is obvious. In Figure 1, the monomer peak is well resolved so the broadening problem is clearest. 

9 The aggregates are progressively less well resolved so they show the problem to a lesser extent. 
10 

11 In order to compensate for these distortions, a variety of methods have been employed. Most are 

12 based upon the work of L. H. Tung presented in his 1969 paper published in the Journal of 

13 Applied Polymer Science, volume 13, p 775 et seq. The afore-referenced book by Yau et al 

14 discusses some of the means to correct band broadening. Throughout the chromatography 

15 literature, many papers address means to perform these corrections, Le. to restore the broadened 

16 band to the shape it would have had were no broadening sources present. Most of the techniques 

17 are numerically unstable and can result in physically unreasonable results, such as negative 

18 concentrations, or negative scattered intensities, and therefore are rarely used. The 

19 characterization of proteins by MALS measurements represents an area of particular importance, 

20 yet even here, band-broadening corrections are rarely seen because of the difficulties associated 

2 1 with their implementation. 
22 
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Band broadening effects occur in a variety of multi-detector implementations of liquid 
chromatography. For example, if a sample is to be measured by an on-line viscometer so that its 
intrinsic viscosity may be determined, then a dRI will also be needed. As the sample passes from 
the dRI to the viscometer, it will experience band broadening resulting in a spreading of the 
specific viscosity curve with respect to the concentration curve derived from the dRI 
measurement. Again, the computational difficulties, associated with restoring the broadened 
signal to conform to what it would have been without broadening, often result in questionable 
results, especially if the samples are indeed monodisperse such as unaggregated proteins. 

It is the purpose of the present invention to provide an analytical method by which band- 
broadening effects may be corrected in software with both ease and precision. Another objective 
of the invention is to provide a new approach to such corrections whose departure from the 
conventional approach is sharply distinguished therefrom. The present invention is expected to 
be of greatest utility for the analysis of separations of protein samples. A consequence of this 
application will be the more precise characterization of protein conjugates as well as aggregate 
states. A further application of this invention is its ability to improve the determination of 
intrinsic viscosity measurements by correcting the band broadening effects associated with the 
sample passage between a concentration detector and a viscometer. An important application of 
the invention will be its application to the rapid determination of the 2 nd virial coefficient of an 
unfractionated sample as discussed, for example, in U. S. Pat. No. 6,41 1,383 and the application 
10/205,637 filed 24 July 2002, now in process of issuance, by Trainoff and Wyatt. The 
successful implementation of that method depends critically upon the ability to determine 
accurately the square of the concentration at each elution volume. This requires an accurate 
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measurement of the concentration at each such elution volume. The latter determinations are 
made generally on unfractionated samples whose detector responses are a single peak whose 
mass composition is the same at each interval collected. This single peak is broadened as it 
proceeds through the series of detectors leading to systematic errors in the derived results, unless 
the effect of broadening is corrected. 

Brief description of the drawings 

Figure 1 shows an uncorrected plot of the molar mass of a BSA sample with the dRI detector 
downstream of MALS detector. 

Figure 2 shows the hybrid broadening model fitted to the 90° light scattering trace of a 5ul 
injection of 200kD polystyrene in toluene. 

Figure 3 shows the data of Figure 1 corrected by the method of the invention. 

Figure 4 shows an uncorrected plot of intrinsic viscosity versus elution volume for a BSA sample 

with a dRI detector downstream of viscometer. 

Figure 5 shows the data of Figure 4 corrected by the method of the invention. Trace 11 is the 
corrected intrinsic viscosity 

Figure 6 shows an uncorrected plot of the molar mass of a BSA sample with a UV detector 
upstream of the MALS . 

Figure 7 shows the data of Figure 6 corrected by the method of the invention. 
Summary of the invention. 

The present invention comprises a method to correct for band broadening effects occurring 
between detectors of a liquid chromatograph. Rather than correct the band broadening produced 
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1 at a second detector by restoring the band to the form of the first, the method broadens the band 

2 of the first detector to conform to the broadened band at the second. This results in a simpler, 

3 more robust, and more easily implemented means to correct band-broadening effects at all 

4 detectors. Although the method can result in reduced peak resolution when peaks overlap 

5 partially, areas of major application are expected to be for monodisperse samples with well- 

6 separated and resolved peaks. The invention is applicable as well to any number of serially 

7 connected detectors. 
8 

9 Detailed description of the invention 

10 

1 1 Molecules in solution are generally characterized by their weight averaged molar mass M w , their 

12 mean square radius (r g 2 } , and the second virial coefficient A, . The latter is a measure of the 

1 3 interaction between the molecules and the solvent. For unfractionated solutions, these properties 

14 may be determined by measuring the manner by which they scatter light using the method 

1 5 described by Bruno Zimm in his seminal 1 948 paper that appeared in the Journal of Chemical 

16 Physics, volume 16, pages 1093 through 1099. More recently, the method discussed in co- 

17 pending Application 10/205,637 by Trainoff and Wyatt filed 24 July 2002, in process of 

1 8 issuance, represents an advanced technique that will replace the more traditional Zimm approach. 
19 

20 The light scattered from a small volume of the solution is measured over a range of angles and 

21 concentrations. The properties derived from the light scattering measurements are related 

22 through the formula developed by Zimm: 

23 R{0) = K*M w cP(0)[l - lA^ijcPm + 0(c 3 ) , (1) 



Trainoff Band-broadening J 



1 where R (0) is the measured excess Rayleigh ratio in the direction 0 per unit solid angle defined 

2 as R (6) = [/, (0) - / solv (0)] r 2 / I 0 V , /, (0) is the intensity of light scattered by the solution per 

3 unit solid angle, 4, lv (0) is the intensity of light scattered from the solvent per unit solid angle, 

4 /„ is the incident intensity, r is the distance from the scattering volume to the detector, V is the 

5 illuminated volume seen by the detectors, P( 0) is the form factor of the scattering molecules 

6 defined as P(0) = KmR(0)/R(O) , K* = An\dnldcfn\ l(NX) , K is Avogadro's number, 

7 dn I dc is the refractive index increment, n 0 is the solvent refractive index, c is the concentration, 

8 and /io is the wavelength of the incident light in vacuum. The incident collimated light beam is 

9 assumed vertically polarized with respect to the plane containing the light scattering detectors. 
1 0 The form factor is related to the mean square radius (r* } by 

U P(ehl- l -^(r;)sm^0/2) + o(sin\0)/2) (2) 

1 2 Multiangle light scattering instruments measure the amount of light scattered from a fluid sample 

13 as a function of both angle and time. For molecules very small compared to the wavelength of 

14 the incident light, light scattering measurements may be restricted to a single scattering angle, 

1 5 such as 90°, since there would be no measurable variation of the scattered intensity with angle, 

1 6 however in this limit, the molar mass may be measured, but the mean square radius may not. 
17 

1 8 When the instrument is used to analyze the eluant of a separation system, such as size exclusion 

1 9 chromatography, SEC, the sample composition and concentration change with time. 

20 Furthermore, SEC columns typically dilute the sample concentration so that the second and third 
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1 terms on the right hand side of Eq. (1) can be neglected compared to the first. Formally this is the 

2 case when 

3 2^cM„D 1. (3) 

4 Typically, for chromatographic separations combined with light scattering instrumentation, the 

5 above condition is assumed unless is known from prior experiment, in which case such value 

6 may be used in Eq. 1 directly. In addition, if the separation is performed properly, at each point 

7 in time the sample is essentially monodisperse in size, and usually in mass. If the sample is also 

8 monodisperse in mass we have the further simplification that 

9 R(0 9 t) = K*Mc(t)P(0 9 t) (4) 

10 Therefore, the light scattering signal should be directly proportional to the concentration signal 

1 1 and the molar mass can be computed by the ratio 

K c{t) 0_> o 

13 A common test for the correct operation of a MALS system is to inject a nearly monodisperse 

14 sample through a separation column, and then measure both the light scattering and the 

15 concentration signals for a monodisperse peak. If normalized, they should overlay exactly at all 

16 elution times corresponding to the constancy expressed by Eq. (5). However, if the concentration 

17 detector is downstream of the MALS detector, two effects will change the shape of the peak: 

18 mixing and diffusion. When this is the case, the peak appears broadened; the center of the peak is 

19 suppressed and the "wings" are enhanced. Therefore, when the ratio in Eq. (5) is formed, it will 

20 not be constant and the molar mass derived near the peak will be systematically overestimated, 

21 while that near the wings will be underestimated. The situation is reversed when the 

22 concentration detector is upstream of the MALS detector. 
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1 

2 It should be pointed out that samples injected into a chromatograph without fractionation, often 

3 described as flow injection, also experience band broadening during their passage from one 

4 detector to another. Although the molar mass distribution at each successive elution volume are 

5 expected to remain constant, band broadening effects will affect various weight average 

6 properties derived from downstream broadening at selected detectors. Thus an injected sample 

7 peak comprised of a polydisperse molar mass distribution will broaden as it passes, for example, 

8 between a light scattering detector and a dRI concentration detector. When the unfractionated 

9 light scattering signals are combined with the broadened concentration detector signal, the 

1 0 derived weight averaged molar mass near the center of the peak will be calculated to be larger 

1 1 while the corresponding values in the wings will be smaller. This result is exactly equivalent to 

12 that seen for fractionated samples. 
13 

14 Sample mixing primarily depends on the geometry of the system: the length and diameter of the 

1 5 tubing, the number of unions, the roughness of the surfaces, and the number of inclusions caused 

16 by mechanical defects, and the geometry of the measurement cells. The effect of mixing is 

17 primarily independent of the sample composition in the limit that the bulk fluid properties, such 

1 8 as viscosity, are not substantively changed by the sample. Diffusion, in contrast, depends 

1 9 explicitly on the sample properties. 
20 

21 It is instructive to estimate the relative importance of diffusion and mixing for a sample flowing 

22 between two detectors. For uniform spheres the Stokes-Einstein relation 

23 A- =7^' (6) 
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1 relates the diffusion constant D T to the sphere radius r, where k B is Boltzmann's constant, T is 

2 the absolute temperature, and 77 is the solvent viscosity. 
3 

4 There are many sources of mixing, but to estimate the scale, consider the simple model of a well- 

5 mixed reservoir. Assume that the light scattering flow cell is a mixing chamber with a volume V. 

6 Let the flow rate be denoted by f, and assume that the concentration of the sample entering the 

7 chamber is c g (t) . The concentration of the sample exiting the flow cell is given by 

8 c'(t) = [ Ci (t)-c(t)]/t f (7) 

9 where t f = V I f is the time to fill the reservoir. Consider the effect of a narrow pulse, containing 

10 a mass m 0 of sample injected into the flow cell so thatc,(0 = m 0 S(t) . The concentration of the 

1 1 sample exiting the cell is then 

12 c(t) = ^exp(-t/t f )0{t), (8) 

V. 

13 where 6(t) = 0 for t < 0 and 6(t) = 1 for t > 0 . To put this in context, assume that the flow cell 

14 has a volume of V = 80/// and is fed by a flow rate of / = 1 .0 ml I min . The exponential time 

15 constant for the broadening is t f = 4.8 sec . If this mixing chamber is fed and drained by capillary 

16 tubing with an internal diameter of 025mm , the initially localized concentration pulse is smeared 

1 7 out over a length of tubing 1 .6w in length. Of course, when the sample flows through the light 

18 scattering cell, it is not "well mixed," but this is a good estimate of the length scale of the 

19 mixing. 
20 
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1 Next, consider the effect of diffusion for a small molecule. The molecules will diffuse, on 

2 average, a distance 



4 where, D T is the molecule's translational diffusion constant and t is time. Again, to put this in 

5 context, consider the protein bovine serum albumen, BSA, which has a diffusion constant of 

6 7.7 x 1 0~ 7 cm 2 1 sec . How far does the concentration spike diffuse in 4.8 sec ? Equation (9) 

7 indicates that it will diffuse along the tubing a mere 27 jum ! Therefore, we conclude that the 

8 dominant mixing mechanism is non-diffusive in origin and therefore should not depend on the 

9 diffusion constant, and hence the size, of the molecule. The dominate mechanism of inter- 

10 detector broadening is due to mixing and is sample-independent. Therefore if the geometrical 

1 1 broadening of the system can be characterized with a reference sample, the resulting broadening 

12 parameters can then used to correct all subsequent data runs. 



14 Consider a simple linear model for interdetector broadening. For the purposes of discussion, 

15 assume that the concentration detector is located downstream of the MALS detector. 

16 Furthermore, assume that Eq. (3) is satisfied, and that data for a monodisperse reference sample 

17 has been collected. Therefore, the measured concentration depends on the underlying 

1 8 concentration in the MALS detector via the convolution with a parameterized broadening 

19 function5(a 0 ,a 1 ,-",a JI ,r). One has 



3 



{\x\) = JlD~t 



(9) 



13 



20 




(10) 
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1 where c m (t) is the measured concentration at the downstream detector following broadening, 

2 and c(t) is the concentration of the sample when it passed through the light scattering detector. 

3 The parameters a n are model specific and include the width of the broadening function and the 

4 interdetector delay volume. The goal of the broadening correction is to find the best-fit values of 

5 these parameters. According to Eq. (5) one may write 

6 ^(0 = ^£^(0 5 ^0 5 K^i^s^ ? r)Jr (11) 

7 where R (0,t) = lim R {0,t) . Here we have used the fact that the peak is monodisperse, so that the 

8 molar mass M does not change across the peak, and it can be brought outside the integral. In 

9 order to determine the best fit parameters, one forms the chi-squared difference between the 

10 downstream measured concentration signal, c m (t\ and the broadened upstream light scattering 

1 1 signal and integrates it across the monodisperse sample peak 

12 z(^ 09 -,a tt rT) 2 = Jf (^(0-^7 £ *(<U - t)B (a 0 , • V <V^} dt. (12) 

13 One may then use a standard nonlinear least squares fitting package to find the best- fit 

14 parameters for K*M and the a { that minimize % 2 • Denote these best fit parameters a[ . 
15 

16 As a practical matter, one may include other system specific parameters in the above fit. For 

1 7 example, the excess Rayleigh ratio in Eq. (12) is defined as the difference between Rayleigh 

1 8 ratio of the solvent with sample, and the Rayleigh ratio of the pure solvent. The Rayleigh ratio of 

19 the pure solvent can be included as an adjustable fit parameter, instead of manually determining 

20 it by setting a solvent baseline. Similarly, if the concentration detector has a slow drift, then the 

21 offset and slope of the drift can similarly be included in the fit. It is important to note that it is 
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1 essential that the molar mass is constant over the peak so that the inner integral in Eq. (12) 



3 function. If the mass were not constant, it would be impossible to determine uniquely the best fit 

4 parameters a- . 
5 

6 Once the optimal fit parameters have been determined, there are two ways in which they can be 

7 used in the subsequent analyses. One can attempt to "narrow" the concentration measurements 

8 by performing a deconvolution of Eq. (10) and attempt thereby, to answer the question: what 

9 would a concentration detector have measured if it were coincident with the light scattering 

10 detector? Alternatively, one can artificially broadening the light scattering results to answer the 

1 1 conceptual question of what would the light scattering results have been if they were performed 

12 downstream coincident with the concentration detector? We will argue that the former method, 

13 which has formed the traditional approach, is numerically unstable and often gives unphysical 

14 results, while the latter method is numerically stable, but at the expense of reducing the 

1 5 resolution of the measurement. Since the broadening is usually a small correction, the loss of 

1 6 resolution is minimal and the latter method is preferable. 



2 £*(0,f-r)*(, 



a 0 , • • • , a n , r ]dr contains only the measured Rayleigh ratio and the broadening 



17 



1 8 Consider the traditional deconvolution method. One may reconstruct c(t) from knowledge of 



19 c m (t) by writing the Fourier Transform 



20 




21 



where 



Trainoff Band-broadening 14 



1 c ">) = 7^£ e " cm (')^ (14) 

2 is the Fourier transform of the measured concentration, and B^a l \--- 9 a ll ' 9 a>jis similarly the 

3 Fourier transform of the broadening kernel evaluated with the fit parameters found above. One 

4 may then compute the molar mass as a function of time by using Eq. (5). The problem with this 

5 procedure lies in computing the ratio c m (co) I B {a{ , • ■ • , a n ' , <yj . Both c m (co) and 

6 B ^a/ , • • • , a n ' , 0)j tend towards zero as co -> ±oo . Therefore, the ratio tends towards 0/0 which is 

7 indeterminate. Because c(co) contains experimental noise, the ratio will have increasing large 

8 fluctuations for large co . Therefore, when one computes c{t) using Eq. (13), the result contains 

9 high frequency noise and unphysical "ringing," including negative concentrations. While the 
10 noise can be filtered out, the unphysical ringing cannot. 

11 

12 The second method, which is the subject of this patent, is to broaden the light scattering peak. 

13 Therefore, we write 

14 M (t) = _J-^Hm ^R{ej-r)B[a;^^a; 9 r)dr . (15) 

1 5 This constitutes a measurement of the molar mass as a function of time, corrected for 

16 interdetector broadening. Note that right hand side of Eq. (15) contains only measurable 

1 7 quantities, along with the broadening function evaluated with the previously determined 

18 parameters. 
19 
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1 Models of the broadening function 

2 The algorithm described above allows one to correct the light scattering results for the effects of 

3 interdetector band broadening given a parameterized model of the broadening kernel. However, 

4 to apply the algorithm, one must choose a model to use. In this section we will describe a 

5 specific model which works well for the MALS hardware, DAWN instruments, built by Wyatt 

6 Technology of Santa Barbara, CA. However, it is important to note that the algorithm applies 

7 equally well to any definite model. 
8 

9 The light scattering cell used in the Wyatt MALS instrumentation is a cylinder 1 .25 mm in 

10 diameter fed, at 90 degrees, by a capillary that is 0.1 mm in diameter, and drained by a capillary 

1 1 that is 0.2 mm in diameter. Typical flow rates are of the order of Imllmin . One can determine 

12 whether the flow is laminar or turbulent by computing the Reynolds number, defined 

13 as Re = vd I v , where v is the average flow velocity, d is a characteristic dimension of the 

1 4 system, and v is the sample viscosity. For the flow in the inlet capillary one has v = 1 27mm / sec , 

15 d - OAmm , and v = 0.89 cm 2 1 sec , so that Re = 14 . Fluid flow typically does not become 

16 unstable until it reaches Reynolds number of a few hundred. Therefore, we conclude that the 

17 flow in the inlet capillary is laminar. Next, consider the jet of fluid as it enters the flow cell. The 

18 flow rate is the same, but the characteristic dimension increases to d =125 mm. Therefore, the 

1 9 Reynolds number increases to Re = 1 78 . Therefore, we conclude that there is some turbulent 

20 mixing near the inlet of the cell, which has a mixing volume of some fraction of the cell volume. 

21 Similarly, there is some mixing in the inlet of a concentration detector flow cell. 
22 
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1 There is also some instrumental broadening which comes from the fact that the system measures 

2 a finite volume of fluid, and the measurement system has some intrinsic filter time constant. 

3 Therefore the model we will consider is an exponential broadening, to represent the mixing in 

4 the flow cell convolved with a Gaussian term to represent the instrumental broadening. The 

5 model for the broadening kernel is 

6 B ((7,w 9 T) = —^e- T2/2<j2 ®-e- T,w 0(r) 9 (16) 

7 where a is the instrumental Gaussian width, and w is the broadening due to mixing in the inlet of 

8 the cell and the convolution is defined as 

9 a(t)®b(t)= ^a(t-t')b(t')dt\ (17) 

10 To test the appropriateness of this model, we performed the following test. We filled a 5jul 

1 1 injection loop with 200kD polystyrene dissolved in toluene and injected it directly into the light 

12 scattering flow cell. The flow rate was 0.5m/ Imin and the sample collection interval was 0.25 

13 sec. We model the data as a delta function, to represent the injection loop, convolved with the 

14 broadening model. To compare this with the 90° light scattering signal, one forms x 2 as 

!5 z(a,t 0 ,or,w) 2 = [J^R(90\t)-aB((7 9 w i t-t 0 j]dt 9 (18) 

16 where t 0 is the time required for the sample to flow from the injector to the flow cell, and a is a 

17 scale factor to convert from dimensionless units to Rayleigh ratio. There are four fit parameters: 

18 a, t 0 , cj, and w. 
19 

20 Figure 2 shows the results of the fit. The light scattering data points 4 are normalized so the peak 

21 has unit amplitude. The resulting fit to the broadening function is 5. This includes only the 
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1 broadening which occurs between the injection loop and the center of the flow cell, where the 

2 measurement takes place. It is clear the broadening function fits the data well, and therefore is 

3 an appropriate model for interdetector broadening. This broadening model will be used in 

4 demonstrations described below. 

5 Applying the Method to an Arbitrary Series of Detectors. 

6 In this section the method will be generalized for two or more arbitrary detectors connected in 

7 series. Define the detectors responses as D i (t) , as the sample flows serially from D x , to D 2 , . . . 

8 to D n . In general, the sample peaks become progressively broader as they elute through the 

9 system. Any analysis on the data(Z) 1 (f),Z) 2 , (*),...,£>„ (t)) , will be compromised by this 

10 progressive broadening. As before, it is numerically more stable to broaden the upstream 

1 1 detectors to be consistent with the detector with the broadest response. Usually this will be the 

12 last detector in the chain. However, if one of the intermediate detectors has a large instrumental 

13 broadening, the results can be referenced to it instead. For the purposes of discussion we will 

14 assume that the upstream detectors signals will be broadened to match D n . 
15 

16 The procedure has two phases. In the first phase, one determines the broadening parameters by 

17 measuring a monodisperse sample at a low enough concentration that the detector signals, in the 

1 8 absence of broadening, would be directly proportional. Practically, one injects a nearly 

19 monodisperse sample through the chromato graph, and selects the monomer subpeak with which 

20 to determine the broadening parameters. For each detector i from 1 to n - 1 , one computes 

21 *MA.*,.a»)- \J^n{t)-fii]iDXt-r)B{ a ^-T)dTpt, (19) 
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1 where %l (A > r ,>#y ) > is the chi squared parameter for the i ih detector. The parameters p. are 

2 proportionality factors which bring the detector signals to the same scale, x. are the inter- 

3 detector delays associated with the time it takes the sample to travel from one detector to the 

4 next, and a i} are they broadening parameters for the 2 th detector. These x 1 are minimized using 

5 a standard nonlinear least squares fit algorithm to determine the best fit parameters, p[ , z[ and 

6 a- . Software capable of such fitting is readily available in various commercial packages. The 

7 procedures and some of the various software approaches have been discussed in detail by Hiebert 

8 in his 1 98 1 paper "An Evaluation of Mathematical Software That Solves Nonlinear Least 

9 Squares Problems" in the ACM Transactions on Mathematical Software, volume 7, pages 1-16. 
10 

11 In the second phase, one applies the broadening parameters determined in the first phase to 

12 correct all subsequent data runs. In this phase, all of the upstream detector signals are broadened 

13 as if they were coincident to the last detector in the chain. Then they can be compared directly to 

14 the last detector D n . The broadened data are 

15 />,*(*) = £A(^-^(<,r-r/)rfr 5 (20) 

16 Then the original analysis proceeds on the data sets i^D x b (t) 9 ... 9 D n _ x b (t),D n (/)) . 

1 7 Demonstration of the algorithm 

1 8 In this section we will demonstrate the algorithm on sample data for the following instrument 

19 combinations: MALS+dRI, Viscometer+dRI, and UV+MALS. Figure 1 shows data from an 

20 100/// injection of BSA into a Shodex OH pack KW protein separation column manufactured by 

21 Showa Denko, Japan. The flow rate was 0.5ml/min . The fractionated sample first passed 
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1 through the light scattering detector and then the dRI detector. The BSA sample consisted 

2 primarily of monomers, but had a substantial content of aggregates which were separated by the 

3 column. The peak between 16 min and 18 min primarily consists of pure monomer. The peak 

4 between 14.4 min and 16 min is pure dimer, and the subsequent peaks are higher aggregates that 

5 are less well separated. The horizontal axis is time in minutes. Trace 1 is the computed molar 

6 mass M w {t ) divided by the monomer molar mass M wQ . Therefore the monomer reads 1 on the 

7 vertical axis; the dimer reads 2, etc. Trace 2 is the 90° light scattering signal and trace 3 is the 

8 refractive index signal. Both the refractive index and the light scattering have been normalized so 

9 that their peaks have the same maximum height. The vertical scale for traces 2 and 3 is arbitrary. 
10 

1 1 The first step in the algorithm is to choose the broadening model. For the data set in Figure 1, the 

12 hybrid model described in Eq.(16) was used. The second step in the algorithm is to choose a 

13 monodisperse peak on which to determine the fit parameters. The monomer peak between 16 

14 min and 18 min was used. When one forms the z 2 in Eq. (19), there are four parameters to be 

1 5 determined by nonlinear least squares fitting. An extra parameter was added to account for a 

16 possible baseline shift between the two data sets. Therefore the % 2 that was minimized was 

17 z 2 {fi 9 r 0J c7,w,x 0 )= \^[dRl(t)-p^ (21) 

18 where, LS(t) is the 90° light scattering signal as a function of time, and dRl(t) , is the dRI data 

19 as a function of time. This model was minimized by using a commercial Marquart nonlinear 

20 least squares package such as described by D. W. Marquardt in his 1963 article in the Journal of 

21 the Society of Industrial and Applied Mathematics, 11 An algorithm for least squares estimation of 

22 nonlinear parameters," volume 11, pages 431 to 441. Note that to perform the broadening 
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1 correction proscribed by Eq.(20), one does not need /? , or * 0 , so they will not be presented in the 

2 results table below, but they were included to insure that the nonlinear minimization worked 

3 correctly. 



Parameter 


value 


& 


0.0636 sec 


W 


3.178 sec 


t 

de- 


21 sec 



4 These results can be interpreted by noting that a □ w , implies that, as expected, the dominant 

5 broadening effect is mixing, not Gaussian diffusion or instrumental broadening. The r 0 

6 parameter simply indicates that the sample took 2 1 seconds to traverse between the two 

7 instruments. In this context Eq. (20) can be rewritten as 

8 LS b (t)= £z5(r-r)^(or',w',r-r 0 ')t/r (22) 

9 Figure 3 shows the broadened light scattering signal 7 computed using Eq. (22) with the fit 

10 parameters found in the first step above. The light scattering analysis proceeds as before using 

1 1 Trace 7 instead of the original light scattering data. The results are shown in Trace 6. There are a 

12 number of features to notice. First is that Trace 7 now overlaps Trace 2 at the monomer peak. 

13 This shows that model is correct, and that the nonlinear minimization worked properly. The 

14 second observation is that the dimer peak, between 14.5 min and 16 min, is now also a plateau as 

15 expected. The trimer peak, between 13.5 min and 14.5 min is also flatter, but since it was not 

16 baseline separated from its neighboring peaks, the results are less clear. The next observation is 

17 that the plateaus are now integral multiples of each other. This is to be expected since the molar 
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1 mass of the dimer should be exactly twice the mass of the monomer. Similarly, the trimer plateau 

2 is now three times that of the monomer. 
3 

4 Figure 4 shows data from an 100/// injection of BSA into a Shodex OH pack KW protein 

5 separation column. The flow rate was 0.5ml/min . The fractionated sample first passed an online 

6 bridge viscometer and then through a dRI detector. The viscometer measures the specific 

7 viscosity which is defined as 

8 %(t) = rj(t)/rj 0 -l, (23) 

9 where ?](t) is the viscosity of the fluid sample combination and ?] Q is the viscosity of the pure 

10 fluid. The dRI detector measures the refractive index of the sample which can be converted to 

1 1 concentration if the sample dn I dc is known. One can then compute the intrinsic viscosity 

12 defined as 

13 ^M = ^%{t)/c(t) : (24) 

14 If the concentration is sufficiently low then the intrinsic viscosity can be approximated as 

15 *lM*ri sp {t)lc(t). (25) 

16 Trace 9 is the specific viscosity, and trace 10 is the dRI signal. Traces 9 and 10 have been 

1 7 normalized to have the same maximum height. The vertical scale for these two traces is arbitrary. 

18 The results of this measurement are shown in Trace 8 which presents the intrinsic viscosity on 

19 the vertical axis in units of ml/g. The intrinsic viscosity should be constant across each sample 

20 peak. The curvature seen is due to interdetector. broadening, but in this case the viscometer 

21 signal is substantially broader than the dRI signal. This is because the viscometer measures the 

22 viscosity over a much larger fluid volume than the dRI. In this example, the dominant 
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broadening is instrumental broadening internal to the viscometer. In fact, the dRI signal is 
broadened by the passage of the sample between the two detectors, which reduces the overall 
difference in peak widths. The method of this invention can be applied as before, but in this case, 
the dRI signal will be broadened to match the viscosity signal. The first step in the algorithm is 
to choose a broadening model. Again the hybrid Gaussian-mixing model was used. The second 
step in the algorithm is to choose a monodisperse peak on which to determine the fit parameters. 
The monomer peak between 15.2 min and 16.2 min was used. As before Eq. (19) can be 
rewritten as 

^ 2 (/?,r 0 ,o-,w,x 0 )= L(^0-^£^ (26) 

and the fit parameters are determined by minimizing j 2 with respect to the adjustable 
parameters. We find the fit parameters 



Parameter 


Value 


a' 


0.01 sec 


W 


9.04 sec 




15.6 sec 



In this case, the parameter w is 9.04 seconds, which arises not so much from inter-detector 
mixing, but from the large internal volume of the viscometer capillaries relative to the volume of 
the dRI flow cell. The Gaussian term is extremely small, which again indicates the diffusion is 
insignificant. The dRI signal may be broadened by rewriting Eq. (20) as 

dm b {t)^ ^dm(t-r)B(^\w\r-r;yr 9 (27) 
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and the analysis can proceed with the broadened dRI signal Figure 5 shows the broadened dRI 
signal in Trace 12 and the corrected intrinsic viscosity signal is shown in Trace 11. The results 
now consist of plateaus for each peak. 

The last example, illustrated in Figure 6, consists of a lOOjul sample of 3.0mg/ml BSA 
measured by a UV detector followed by a MALS detector. In this case, the LS signal shown at 
14 is significantly broader than the UV signal at 15, giving rise to the curved molar mass trace 
shown at 13. For this example, the UV signal will be broadened to match the LS signal. Equation 
(19) is rewritten as 

Z 2 (0>T O9 a 9 w,x o )= j^LS(t)-fi (28) 

Again the hybrid mixing and Gaussian model is used, and the monomer peak between 13.5 min 
and 15 min is used to determine the fit parameters. Minimizing with respect to the fit 
parameters yields 



Parameter 


Value 


a' 


0.22 sec 


W 


3.61 sec 


t 

*0 


39.5 sec 



Again, the dominant source of broadening is internal mixing. The UV data is now broadened, as 
shown in Figure 7 by computing Eq. (20) rewritten as 

UV b (t)= ^UV(t-T)B(aW,T-T 0 'yT. (29) 
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1 The broadened UV signal is shown in Trace 17, and the molar mass is shown in Trace 16. As 

2 with the previous two cases, the data consists of plateaus, and the effects of inter-detector 

3 broadening are dramatically reduced. 
4 

5 The foregoing demonstration of the inventive method illustrates but a few of the possible types 

6 of band broadening corrections that may be addressed and corrected by the inventive methods 

7 disclosed in this specification. As will be obvious to those skilled in the chromatographic arts, 

8 until this invention, the effects of band broadening have represented significant impediments to 

9 improving the precision of measurements associated with the use of multiple detectors. This 

10 inventive method and all its variants fulfill the expectation that band broadening effects will 

1 1 thereby be reduced dramatically. Accordingly, I claim 
12 

13 
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